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Prafaoa to Third Yearly Report xmder Contract HAS2-7613 


Hopk under Contract NAS2-7613 started on July !• 1973. The 
research goals for a 3 year period stated in tiiis contract are: 

(a) Assess analytically the effects of fuselage motions on st^ility 
and random re^onse. The problem is to develop an adequate but 
not overly complex fli^t dynamics analytical DK>del and to study 
the effects of structin>al and electronic feedback, particularly 
for hingeless rotors. 

(b) Study by conputer and hardware experiments the feasibility of 
adequate pertuzbation models from non-linear trim conditions. 

The problem is to extract an adequate linear pertiubation model 
for the pinpose of stability and random motion studies. The 
extraction is to be performed on the basis of transient responses 
obtained either by computed time histories or by model tests. 

(c) Extend the experin^ntal methods to assess rotor wake-blade 
interactions by using a 4-bladed rotor model with the capability 
of progressing and regressing blade pitch excitation (cyclic 
pitch stirring), by using a 4-bladed rotor model with hub tilt 
stirring, and by testing rotor models in sinusoidal up or side flow. 
Seven reports oi» the work under Contract NAS2-7613 have been 

submitted, references 1 to 7. Reference 1 completes research goal (a). 
References 2, 4, 6, 7 pertain to research goal (b). It is incomplete 
to date. References 3 and 5 pertain to research goal (c). It also is 
Incomplete to date. Reference 3 presents the results of extensive 
frequency response tests. Reference 5 presents dynamic downwash 
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n^asurenants in bovaring during harmonic rotor oxcltatlon. Three 
manuscripts for publication have been prepared. The first hse been 
published as reference 8. Tuo further manuscripts covering part of 
the material In references 3 and 7 have been sxibmltted to journals. 

The extensive rotor state and parameter Identification work with 
computer simulated transients was found to be very tiseful for the 
siibsequent data processing of the measured transients. It allowed to 
sort out possible inadequacies of the spiled identification algorithm 
and of the inputs from possible inadequacies of the measurements and 
of the applied mathematical rotor model. Rotor state and parameter 
identifications from measiired transients are presently complete for 
hovering conditions using cyclic pitch stirring transients. They are 
presented in reference 7. The corresponding work for forward flight 
conditions and for hub tilt stirring is planned to be completed in 
FY 1977 during an authorized extension of the research contract. It 
is also planned to refine the analytical rotor model used for the 
state and parameter identifications to include blade flexibility. 

Rotor state and parameter identification from transients is still 
a field where little experience is available. As elaborated in 
Chapter 1 of reference 4 it takes four important ingredients to perform 
a successful state and parameter identification from transients; a 
suitable input, a suitable instr\imentation for measuring key state 
variables, an adequate mathematical model of the system, and an 
efficient criterion function for the estimation algorithm. In all 
four respects considerable work had to be done in order to finally 
establish a combination of these four ingredients that led to success. 
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Hie leesone leamsd In this effort are ^iped to pave the wa^ for a 
wider ^tlication of rotor dynamic perturbation state and paraioeter 
identifications fr^ traiwiaits about non-linear trim conditions, 
both in rotor wind tunnel testing and in rotor fli^t testing. 
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ROTOR DYNAMIC STATE AND PARAMETER IDENTIFICATION 
FROM HOVERING TRANSIENTS 


Part II of Third Yearly Report under Contract NAS2-7613 


Abstract 

State and parameter Identifications based on a fozm of the maximum 
likelihood method are applied to the problem of extracting linear per- 
turbation models including rotor dynamic inflow effects from transient 
blade flapping measurements. The estimation method is first studied 
in con^uter simulations and then applied to cyclic pitch stirring 
transients generated with a four-bladed rotor model operating in 
hovering trim conditions. The analytical pertvu?bation models extracted 
from the transient test results are compared with transient and 
frequency response tests not used in the state and parameter iden- 
tification. The identified analytical pertm<bation model is also 
compared with a simple theozy. The method that is applicable bolii to 
small scale and full sccile dynamic rotor testing is being extended to 
pertuxt>ations from forward flight trim conditions. 
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aerodynamic blade constant « proportional to y 

Fourier coefficient 

Induced flow gain 

information or sensitivity matrix 

non-dimensional natiuc>al blade fl«^plng frequency 

rotor radius 

blade section lift slope 

rotor-to-ground plate distance 

distance from rotor center 

non-dimensional time (time of revolution 2ir) 

flapping angle, positive up 

blade Lock number 

up blade pitch angle 

down induced flow parameter, or innovation 
standard deviation, or rotor solidity ratio 
standard deviation divided by mean value 
induced flow time constant 
non-dimensional pitch stirring frequency 


Superscripts 

* time derivative 

* estimate 

trim or bias value 

* equivalent value 
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multiblade variables 
mas value 
neaexarad variable 
first Fourier coefficient 
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Introduction 

Chapter 2 of reference 4 gives the results of state and paraineter 
identifications based on computer simulated measurements of transients 
fi^ hovering trim conditions. These simulation studies have been 
improved with respect to input design and data handling. The main 
progress, however, was made in using transient test data rather than 
simulated data for the identification of perturbation models for 
rotor hovering trim conditions. In previous e]q>erimental wovk the 
method of frequency response testing was used to obtain data on 
dynamic rotor inflow, see reference 3. During cyclic pitch stirring 
extensive zero advance ratio dynamic rotor wake measurements were 
obtained at a distance .2R below the rotor plane and are presented in 
reference 5, that also contains a description of the test set-up, the 
planned test procedures and the planned data handling systems for the 
subsequent cyclic pitch transient testing of the model rotor. 

In the present report we discuss improvements in the test equipment 
and data handling made during the report year, present transient test 
results for hovering trim conditions, and compare the identified per- 
turbation model with a simple theory eind with the results of the 
preceding frequency response tests. Some transient test data have 
been obtained in the wind tunnel for forward fli^t trim conditions. 

Data processing for these tests is not conq>leted as yet. 
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Furthw Simulation Studies 

In Ouster 2 of rafirsnoe 4 th« maximum likelihood method of 
state and parameter Identlfioatlon is ei^lained and applied to a 
perturbation model for hovering trim conditions developed in ref> 
erence 8. A step input of cyclic pitdi was applied at time t » 0. 

After the transient from the step input had almost subsided, pitch 
stirring with a constant stirring acoeleratlon was initiated at 
t - 70, and the simulated test data between t s 70 and t « 82 were 
used for the state and parameter identification. It was found that 
better results are obtained for progressing cyclic pitch stirring 
than for regressing stirring. It was also found that the inclusion 
of the Initial values for the fl^ping and Inflow variables in^roved 
the results thou^ at an increase in con^uter CPU time per 
iteration from 3.8 seconds to 6.3 seconds. 

The system and measurement equations for perturbations from 
hovering trim conditions are given in zaference 4 and will not be 
repeated here. We fixat ^plied these equations to the method of 
optimum data utilization developed in Chapter 3 of reference 4. Fig. 1 
shows the relative standard deviations for the 3 parameters A, AL/t, 

1/t vs. time of the transient for a progressing pitch stirring 
acceleration of m = >.l/ir. The true values of the parameters are 
A = .50, AL/t s .25, 1/t s .125. The assumed measurement noise in 
Gaussian with zero mean and a standard deviation of c($g,) = .05. The 
fourth identified parameter P (true value 1.2) was very rapidly and 
accurately identified and is not shown in Fig. 1. It is seen that a 
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translant tlma of t « X2 as used in raferenoa 4 is inadequate while at 
t s 18 the standard deviations for all parameters have reached approxi- 
mately their minimum asyn^totic values. More data will not improve the 
accuracy of the paranmters aiqr s^preciable amount. Spot checks of 
identifications with various data lengths confizmied the characteristics 
shown in Fig. 1. 

In order to avoid the Identification of the initial conditions as 
was done in reference 4, the input was changed such that at the time of 
initiating the cyclic pitch stirring transient the cyclic pitch is zero. 
The input for oi = -.l/ir is represented by Fig. 2 in this report rather 
than by Fig. 2 in Chapter 2 of reference 4. The 6 j input is the same 
except for the extended time of the transient from t = 12 to 18. The 
0JJ input curve is shifted downward by 1.5 degrees in order to begin 
with the value zero. and djj were polluted by computer generated 
Gaussian noise with standard deviation O 0 = .10. The flapping response 
of blades 1 and 2 polluted by conputer generated noise of standard 
deviation ~ .05 is shown in Fig. 3 This figure can be compared to 
Fig. 5 in Chapter 2 of reference 4. The two responses are quite 
different. 

Table 1 gives for a time interval At = .1 used in the numerical 
integrations with the Runge-Kutta method and for a data length of t = 0 
to 18 the results of 3 iterations. The value P = 1.20 weis assumed given. 
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Table 1( P s 1.20 


Parameter 

A 

AL/t 

1/T 

True Value 

.500 

.250 

.125 

Initial Estimate 

.400 

.200 

.083 

Iteration 1 

.487 

.243 

.122 

2 

.492 

.244 

.121 

3 

.492 

.244 

.121 

-1/2 

M ' /Ident. Param. 

.028 

.13 

.14 


The iterations converge very rapidly to parameter values that are only 
little different from the assumed true values, taken from reference 8. 
Rapid convergence and small parameter erixirs were also found when the 
initial estimates differed much more from the true values than in the 
preceding table. The CPU time per iteration on the IBM 360/55 
computer is 3.8 seconds, down from 6.3 seconds in cases with initial 
value identification in Chapter 2 of reference 4. Because of the 
modified input shown in Fig. 2 initial vedues of the perturbation 
state variables Bj, Sjj* Vj* are zero and need not be identified. 

The relative standard deviations for the parameters indicated in the 
last row are somewhat larger than in Fig. 1 for t = 18. The reason is 
that now 6j and were also noise polluted, see Fig. 2, while Fig. 1 

assumed noise pollution only in the simulated flapping measurements. 

It also is of interest to identic the equivalent Lock number. 

Using the inputs of Fig. 2 and the simulated measurements of Fig. 3 
one obtains the following 3 iterations 
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ParoMter 

A* 

Initial Estimate 

.400 

Itex^ion 1 

.280 

2 

.272 

3 

.271 


Thus the equivalent Lock maaber is .271/. 5 = .54 times the true Lode 
nuid>er. 

Fig. 4 shoifs a comparison of the flapping responses for zero 
dynamic rotor inflow, for a ^rnamic inflow represented by the 
equivalent Lock number of .54 times the true value, and for the 
correct rotor dynamic inflow representation by L and t. (A = .500, 
AL/t = .25, 1/t = .125). Using the equivalent Lock nuidt>er does 
improve the response as conpared to that without inflow. However it 
does not properly represent the dynamic rotor inflow effects. Finally 
Fig. 5 shows the responses of the cyclic inflow variables v^, to 

the input of Fig. 2. These values cure quite large and explain why 
omitting dynamic inflow leads to large response errors. 

Improvements in Test Equipment and Data Handling 

The acceleration for the <qrclic pitch stirring electric motor had 
originally been measured in the horizontal position of the shaft. 

When the first series of 30 tests at 2°, 5®, and 8® collective pitch 
setting was evaluated it was found that the maximum likelihood identi- 
fication method did converge as it did in the simulations. However, 
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the obtained parameters were in error. The blade Lock nuedser converged 
to approximately its equivalent value rather than to its acttial value* 
idiich made the values for L and t meaningless. It was then 
discovered that the pitdi stirring acceleration used for the tests 
was much smaller than anticipated from the preceding measurements. 
Apparently* changing the pitch stirring motor axis to the vertical 
position used for the pitch stirring tests largely reduced the 
acceleration. 

Conq>uter simulations were then conducted for the low pitch 
stirring acceleration of the test (m = .05/ir) and a duration of the 
transient from t = 0 to t = 12. The same identification difficulties 
occurred as for the test data* proving that the cause had been an 
inadequate input. Later it was found that low pitch stirring 
accelerations are adequate, if the length of the trcinsient used in 
the parameter identification was sufficiently lengthened. The pitch 
stirring motor was then remountf ’ so that the pitch stirring shaft is 
driven by a timing belt drive. The drive has a transmission ratio up 
to 1:7. Fig. 6 shows the input obtained for two pitch stirring 

accelerations iised for the transient tests to be described in the next 
section. The identification length was extended to about 23 time units* 
since the pitch stirring acceleration is non-uniform and lower than 
average at the beginning of the transient. 

Simultaneously with the modified pitch stii?ring drive the 
change-over from two differential blade flapping meeisurements to four 
separate blade flapping measurements was made. This required the 
installation of a new 20- ring slip ring assembly. As mentioned in 
reference 5 this modification had been planned for the forward flight 


transient tests 
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Data pvooaaaiag of tiis tvanaienta ravaalad iaaeciiraoioa in tiie 
aaiimitii oaaatireiiBntB of pitch stirring and rotor shafts. Thasa 
Inacoureoias nara nodnly causad by the rssolsar wavs shtgias that 
svhstantially dariatad from a haroxinic function. After rsshsping 
the rasolrsr wavs form by ei^irically finding an optimum ftltar and 
excitation frequency the azimuth errors ware reduced. Aamoving an 
electric interference between the resolvwe fUz'ttier iiqiroved tl»ir 
accuracy. 

A rather strong 2 per rev. oscillation was observed in the 
hovering blade flapping measurements prior to the transient. It was 
found after cox^i^rable ei^rimentation that this phenoaenon <x>uld be 
8td>8tantially reduced by correcting non->uniformities of pre-lag and 
pre-fl^ angles between the individual blades. 

As compared to the procedme described in reference S a number of 
DKMlifications were made in the data processing. The entire data 
processing and state and parameter identification program was 
reformulated for use on the PDP-8 and POP-12 mini coiq>uter complex 
that was programmed to perform all necessaxy data manipulations; 
digitizing the analogue data» transformation from rotating to fixed 
coordinates* separating the transient pertwbation values from the 
trim values, executixig the maximum likelihood algorithm to (htain the 
unknown parameter estimates, and generating anal]rtical responses based 
on the identifiev' parameters. This avoids manual key punching and is 
also more cost effective, since 20% usage time for the PDP-8 and PDP-12 
mini computer complex could be bou^t under favorable conditions. The 
core nemory of the complex had recently been increased from 16k to 
38k which proved to be adequate for tiie task at hand. 



10 


The idsntificatlon progron was modified to to Hie 

peramstere A, L, t separately* rather than in oonbination. 
Although the individual iteration requires someidiat Icmgsr ocoiputer 
time* the convergenos was found to i^rove* thus saving overall 
conqiuter time. The identification program was furtiier modified tqr 
considering the imiltiblads fla^qilng vnpiables as meraured gjuantities 
and Iqr introducing unknown biases to be identified for these 
variables. Thus the innovation is now given by 


V 



®IIm 




( 1 ) 


where 8^ and are the unknown biases. Ibis method has the 
advantage that the transfbraation from single blade flapping angles 
to raultiblade variables needs to be performed only once for the 
measured quantities, and not for each iteration. 

The identification algorithm was expanded to include the parameter 
P(blade flapping frequency when rotating). In Ihis form all parameters 
in the analytical model are identified. It was found that such a 
coiq>lete identification from test transients is possible and leads to 
P> values within a few percent of the predicted values. For the 
routine identification runs, P was then assumed to be given. 

The identification algorithm was further modified to include for 
each iteration the updated measurement error covariance matrix obtained 
from the differences between the experimental data and those predicted 
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wi^ the psranstttVB for the prooediag Itaratlcmt It was found that 
batter eonverganos is obtained if only the diagonal terns of the 
neasurensnt enrar covariaaoe are used as "aai^ts" for 
criterion function used in the subsequent iteration. Ihe i^dated 
msasurensnt error covarianos matrix ia used also in the determination 
of the Cramer-Rao lower bou^ H for the parameter estimatest 
which thus becomes mors meaningful. 

The effect of transient duration tiie parameter id^itificatlons was 
not only studied for the simulatioiffi but also for the tests. It was 
found that for good convergence and data fit a minimum data length of 
about 23 time units is required as indicated in Fig. 6. This is 
scmiewhat more than the 18 time units necessaiy according to Fig. 1 
for simulatioQs for a c 7 clic pitch stirring acceleration of u s -.1/x. 

Transient Test Results 

The test conditions extended over the entire collective pit<h 

range iqi to stall and covered three TOtor-to-groimd plate distances of 

1.28R, 1.02R» and .78R. Only one exanple of intermediate steps in 

the identification process will be given here in form of Table 2. The 

collective pitch setting is 6.3°, the time interval between data points 

is .1. As indicated in Fig. 6 the identification was performed with 

the transient between t = 0 and t = 23. The blade flapping frequency 

is P = 1.17. The last column ® measU3?e of 

the overall agreement between analjrtical and test data. The last row 

- 1/2 

gives the diagonal terms of M indicating the parameter standard 
deviation (CrameV-Rao lower boutnd) for iterations 2 to 4. 
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Iteration 

L 

Td>l« 2 

T 

P » 

m 

1.17, 

*n 

®0 * 
A 

6.34 

0 

6.00 

8.00 

0 

0 

.600 

.0668 

1 

4.30 

3.57 

.025 

.025 

.395 

.0276 

2 

4.35 

4.38 

;026 

.023 

.394 

.0253 

3 

4.46 

4.59 

.026 

.022 

.392 

.0253 

4 

4.44 

4.60 

^2&. 

.022 

.391 

.0253 


.08 

.16 

.001 

.002 

.007 



The accuracy of the anal 3 rtical as cospared to the test data, 

((o|^ + , does not change beyond the second iteration, 

though small parameter changes still occur. 

Fig. 6 shows the input whereby the solid line refers to the 

case of Table 2, the dash line to a hi^er cyclic pitch stirring 
acceleration. The 6^^ values of Fig. 6 are normalized to an amplitude 
of one. The actual cyclic pitch amplitude was 1.4®. Fig, 7 shows for 
the time from t = 0 to t = 23 the flailing responses and ^11 to the 
normalized input from Fig. 6. The conplete analytical model with 
identified A, L, t (solid line) fits the test insults (crosses) very 
well. An Identification of the equivalent A = A* resulted in the 
value A* = .19. The 3?esponse for this simplified anal 3 rtical model 
(dash-dot line) shows some eirors as conpared to test data. Merely 
omitting the dynamic rotor inflow amd retaining the correct value of 
A leads to a flapping amplitude (dash line) that is about twice the 
correct value . 


kJ 
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In addition to ofMopasittg the analsrtical and test vesp«aeea fen the 
trsnatent uiod in the panmaeten identifioation, it ie of inteioat to 
pevfMnn aiuh a oonpariatm also for otiier treBsimte. Fig. 9 refers to 
a trensiMit witii a higher aeoeXeration for which the iiq>ut ia ^own in 
Fig. 6 dash line. Mhile the trend is the same as diseussed tor 
Fig. 7, the agreement between analytical and test data is less 
satisfactoiy. In part ‘diis is caused by the fact that tiie bias 
ezvors for the msasurements have not been reoraved as they were for 
the preceding test. The discrepancies become greater as the qrclic 
pitch stirring frequency increases. This may be catmed by the 
neglected unsteady aerodynamics. The transient shown in Fig. 8 
extends to a ^clic pitch stirring frequency of about 3 times the 
rotor rotational frequency in the fixed coordinate system. It is 
d>viou8ly not possible to cover this entire frequency range with the 
same set of L and t parameter. This becomes evident also From Fig. 9 
that shows the ccanparlscm between anal 3 rtlcal and test results for the 
slow transient bey<md the time 23. The test data in this time region 
have not been used for the paran»ter identification and the ftpequencles 
in this region ar<e much hi^er than in the time region from 0 to 23 
units. Fig. 10 shows the computed induced inflow and for the 
slow transient. A comparison with test data is not possible, but this 
figure gives an indication of the substantial magnitude of the (tynamlc 
Inflow that is of the same order of magnitude as the ii^uts 6^ and 6^^. 

In all transients the parameter A was identified close to A « .39. 
The results for L and t are shown in Fig. 11, lhat summarizes all test 
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rMults obtained. The af^ct of the rot or- to- ground plate distance is 
masked by the scatter in tile identified 1 and t values. 

Compariseo with Theory 

A simple momentum tiieovy of ^^niamic rotor induced flow is pre- 
sented in rnferenoe 9. It is based on the assumption tbat the axial 
induced flow is unifcnnn over the disk and that the harmcmic induced 
flow varies linearly with the radial distance &<om tiie rotor center. 

The air mass participating in the angular pitches and rolling motion is 
determined A?om potential flow theozy about a solid circular disk. 

Under these assumptions and under the further assunption of snmll 
linear perturbations about trim one obtains for sero advance ratio 
the following ei^resslons for L and t: 

L = aa/(2Y 7^) , t * .IWvq , . (2) 

ao/2Y = .113. Thus, according to the theory* 

L = T = . 113/ Vq (3) 

The average time mean Induced flow was measured for 3 collective pitch 
angles in a plane .12R below the rotor, using hot wire anemometzy. 

The results and the L and t values from Eq. (3) are given in the 
following Table 3. 


Table 3 



Vo 

II 

2 

.016 

7.0 

5 

.022 

5.1 

8 

.033 

3.4 


*The relation between 7 q and from the momentum theory is not satisfied 
and needs corrections from gzv?und effects and from 3-dlmensional flow 
effects (blade tip losses). 
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By oonf>arl&g thasa values with those of Fig. 11 It is seen that for 
Qq * 8® there is good agreement between the theoretloal values of L 
and T and those identified from transients. For 0^ ■ 2® and 5® the 
theoretical values for t are too high. The equality of L and t 
demanded by the theory also is not valid over the entire collective 
pitch angle range. Nevertheless , the single theory can be used to 
establish approximate values for the rotor induced dynamic flow 
parameters. For example with L = x = 4.5 for 9^ « 5° one finds that 
the eigenvalue for the regressing flapping mods is changed by the 
dynamic inflow from -.20 ^ 16j to -.06 .14j, while the eigenvalue 

for the progressing flapping mods is hardly a^ected by the ^fnamic 
induced flow. Thus the damping of the regressing flapping mode is 
reduced by two thirds. Such changes are not limited to hingeless 
rotors with hig^ P values and should have a noticeable effect on 
precision hovering. Reference 10 gives more data on eigenvalues. 

Coim>ariaon with Frequency Response Testa 

Extended frequency response tests with constant pitch stirring 
frequencies were conducted with the fouu^b laded rotor model, see 
reference 5. Figs. 12 a to c show for O© = 2®, 5®, 8® respectively 
comparisons between analytical results based on the identified 
parameters and fi!<equen<y response test results. The ateolute values 
of the first harmonic Fourier coefficients for the Induced flow 
(upper portion of Figs. 12 a, b, c) and for the flapping response 
(lower portion of Figs. 12 a, b, c) are plotted vs. Ihe cyclic pltdi 
stirring frequency, -hd values correspond to regressing, -u to 
progressing cyclic pitch stirring. 
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An example of dynanio induced floM meaeuMinente nade at 7 tu^al 
etatiena .12R belox the rotor including the soatterhande (j|^ 1 
standard deviation) is in Fig. 13. Ihe assunption of a linear 

radial distribution is not itell satisfied. The strai^t 1^ 
distributi<Hi ^e«nM in Fig. 13 has l^e same nomsnt about tiie ro£or 
center as the measured distribution. Su<di strai^t lines now UMtd to 
generate the data points (crosses) in the i^per portione of Figs. 12 a» 
b* c. Tbou^ no flow msasurements were used to establish the solid 
curves in Fig.s 12 a* b, c and thou^ the presented flow data (crosses) 
are derived A?om inaccurate measurements* see Fig. 13* the agreement 
in the dynamic flow Fourier coe^icients is reasonable. The 
ei^eriraental A?equen^ flapping responses shown in the lower portions 
of Figs. 12 a* b* c agree even better wilb the analytical results. 

Future Work 

Research Contract NAS2-7613 has been extended through FT 1977. 

Ihe following items* failing under either Research Goal (b) or 
Research Goal (c) listed in the Pzeface* are planned to be studied. 

1. Cyclic Pitdi Stirring Rang) Inputs. A mechanism has been built to 
impose qrclic pitch stirring ramp Inputs. After Idie release of a 
spring the qrclic pitch stirring shaft is turned by an angle of 
90 to 330 degrees. An energy absorber prevents the rebound of 
the shaft after impact on a st<^. It is intended to conpare the 
measured blade fl^ping response for this ranp input with the 
response conputed with the identified pez^uibation model. Several 
hovering trim conditions and possibly forward fli^t trim conditions 
will be used. 


hi<3producibility of i - 

Sal page is 



17 


2. State and Parameter Idwatlfloatlons flpom Forward Fllriit Trim Conditions « 
Oontraxy to hovering conditions • eyolio and mllectiva flying are 
cooked in forward fUi£ht. Hie question is whether tiie eoiq^llag 

is strong enough at low rotor advance ratio to identify states 
and paramstere in tiie collsotivs flawing pertuz4>atiMi equation 
from cyclic pitch stirring transients. Collective pitdi transients 
may be needed, which are not feasible! with liie present wind 
tunnel rotor model. At higher advance ratio the simulations have 
shown that adequate collective - tgrclic flapping coupling is 
available to perform a complete parameter Identification from 
(qrcllc pltdi stirring transioats. 

3. Extend Transient Testing to Hi^er Advance Ratio. 

At and beyond an advance ratio of .7 the tip palii plane was 
observed to show indications of an Intending instability. A new 
excentric with .5° amplitude has been made that can be substituted 
for the present 1.5° amplitude excentric. It is planned to 
retest hi^ advance ratio conditions with the smaller excentric. 

4. Extend the Perturbation Model to High Advance Ratio. 

The present perturbation model to be used for the state and 
parameter identifications neglects effects of reversed flow and of 
blade stall. In the hi^ advance ratio region at least reversed 
flow effects must be included. It is likely that at high advance 
ratio the time constants for djmamic inflow in pitch and in roll 
may be different. The identification algorithm will be extended to 
the case of two separate cyclic inflow time constants. 
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5. Cofflplate CoBctruotion of Hvab Stiratag Modal* TSm oonstvttetioD of 
a hxb ativriag model is almost ooaiploted* It mads to be obeoked 
out askd oalibmated* 

‘foot Hob Stirring Model. Tbs hub stlxving model will produce 
different blade flying transients witii hi^r ooii^ling between 
mlleetive and oarclio fleeing. It is planned to test the model 
both for hovering end for fomard fliid^t trim conditions. 

7. Elastic Blade Perturbation Modal. So far Ibe blades in '^e 
perturbation model have been assuonsd to be rigid and flexibljr 
binged at the rotor center. The blade natural flying ftrequenqr 
could be easily and accurately IdentiHed firom transient tests. 

It is planned to study at least by simulations some means to 
extend the identification method to elastic bla^s. 

8. Sinvffloidal Wind Tunnel Flow. A wind tunnel test section has been 
built that allows to produce sinusoidal i^-or-'down flow by 
oscillating vanes. A preliminaxy flow survey has been performed. 
If time and remaining funds allow, this tunnel test section will 
be used to produce transient normal flow variations at the rotor 
model. These transients could be used to IdentiQr the parameters 
in the collective flapping equation that are difficult to obtain 
from <^clic stirring transients at low advance ratio. 
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Coaolwlons 

1* The state and parameter Ideatifleattm netiied prevleualy seleeted 
(maximum likelihood Identlfloatlon) on the basis of oMoputer 
simulated measursments oould be sinpUfled bp using pertuifbatioas 
from trim with sero initial eanditions of the states. Previously 
the initial ^olio pitdh had been oonsidsred a pertuebation from 
trim, which required Initial value identifiLoations. 

2. An analytical model d^th identified equivalent Look nunber shows 
much better agreement of the tz'snsient responses with the actual 
responses than an analytical model «nitting rotor ^rnando Inflow. 
However* considerable response errors remain. 

3. Both the test equlpmnt and the data processing routine needed 
improvements. Non-uniformities in precone and prelag between tiie 
blades had to be removed. The resolver waves had to be reshaped. 
The cyclic pitch stirring acceleration c£^abiUty had to be 
Increased. Data handling could be Inproved by basing the 
identification algorithm on multiblade data rather than on single 
blade data. Measurement bias for the multiblade variables was 
additionally identified. The entire data processing and identi- 
fication program was reformulated for use on the PDF- 8 and PDF- 12 
mini computer con^lex with increased core memory, so that manual 
card punching was eliminated. 

4. Minimum data length for the transients was found to be a critical 
requirement not only for tiie conputer simulated but also for the 


actual measurements. 



5. After <wr(9riately modifying tiie test equipment, tiie input, and 
the data banding programs, the state and parameter identifications 
from oyolic pitch stirring trmisimts fipom hovering trim oonditiams 
were highly euocessfCl. Convergence was rapid, the agreement 
between predicted and actual response very good, except fm cases 
where Ihe transients included high frequencies whidi the simple 
mal3rtiofid model that is based on quasi-steady aero^mamios 
cannot properly represent. 

6. Comparisons of the identified analytical perturbation models 
ftom hovering trim conditions with a simple theory showed 
reasonable agreement. 

7. Con^ariscxiB cf the identified perturbation models with frequency 
response test results showed good agreement in the flapping 
responses and reasonable agreement even in the quantitative 
description of the dynamic rotor inflow characteristics. 

8. Future work will be mainly directed toward forward flight 
transient cgrclic pitch stirring tests, toward hi^er advance 
ratio conditions, toward hub stirring tests, toward the inclusion 
of blade elasticity, and possibly toward the testing of the rotor 
model in wind tunnel flow transients. 
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n«ii» CaptlqoB 

fig. 1 Iformaliaed Standard teviations va. Oaretim of Qyolic Piteb 
Stizvlag Tranaientt u • 

Og ■ .05» A « .SO, AL/t a .SS. 1/t s .125 

Fig. 2 Noise PoUated Simulated (^roUo Pitsli Stindng witSi 

Angular Acceleration & « -.10/ir, s .i 

Fig. 3 Noise Polluted Simulated Blade Flapping Heasurements for 
m * -.10/», Op * .05 

Fig. 4 Cyclic Fl^^ing ResponsM to Inputs fi?om Fig. 2 

Fig. 5 Cyclic inflow Responses to Ii^uts from Fig. 2 

Fig. 6 Test Ii^ut fOT Two C'clic Pitdk Stirring Acceleratiois 

Fig. 7 Response for Slow Cyclic Pitd) Stirring Acceleration. 

Solid Line: A « .39, L s 4.4, T > 4.6; 

Oesh-dot-line: A » A* « .10, L * t » 0; 

Dash line: A ~ .39, L s t b O; 

Crosses: test points 

Fig. 8 Responses for Fast Qrclic Pitch Stirring Acceleration. 

Lins 83 rnhol 8 same as in Fig. 7 

Fig. 9 Resp<»i8es for Slow Cyclic Pitdi Stirring Acceleration 
Beyond t » 23. Line symbols same as in Fig. 7 

Fig. 10 Induced Flow Variables ft>om Identified Rotor Model with 
A B .39, L B 4.4, T = 4.6 

Fig. 11 Identified Induced Flow ParEDDSters vs. Collective Pitch 
Angle 

Fig. 12 a to c 

First Harmonic Fourier Coefficient Amplitudes for Induced 
Flow and for Flapping Responses, 6^ = 2<*, 5°, 6° 
respectively, P = 1.17 

Fig. 13 Radial Distribution of First Harmonic Induced Flor for 
b) = -.067, - 5® at .12R below Rotor Plane 
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Fig. 13 
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/^pendlx A 

List of Purdiased and Borrowed Equipment 


I Equipment purchased under Contract NAS2-4151 and carried over to 
Contract NAS2-7613 


1. Items costing $500. to $1000.- 

1 TSI Model 1057 Signal Conditioner 
1 TSI Model 1125 Calibrator 

2. Items costing $1000.- and over 

1 TSI Thermo Systems Inc. Model 1050 
Research Anemometer 

2 TSI Model 1052 Polynomial Linearizer at 
$1,150.- - each 

i TSI Model 1015C Correlator 


II Equipment purchased under Contract NAS2-7613 

1. Items costing $500.- to $1000.- 

1 EMI Laboratories 1540-9600 
(Part of Stroboscope) 

2. Items costing $1000.- and over 

None 

III Borrowed equipment over $1000.- New Value 

1 TSI Single Channel Hot Wire Ananometry System 

2 Tectronix Type 502A Dual Beam Oscilloscopes 
1 Ampex FR-1300 Recorder/Reproducer 

1 Hewlett Packard Thermal Recording System 
1 Consolidated Electrodynamic.- Corp. 

Type 5-124 Recording Oscillograph 
1 2-Channel Adjustable High-Low Pass Filter 


$ 645.- 

575.- 


1,285.- 

2,300.- 

1.075.- 

$5,880. 


565.- 


3,600,- 

2 , 200 .- 

7,500.- 

7,500.- 

3,000.- 

2 , 200 .- 

$26,000. 



